Plasma formation from ultracold Rydberg gases 
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Recent experiments have demonstrated the spontaneous evolution of a gas of ultracold Rydberg 
atoms into an expanding ultracold plasma, as well as the reverse process of plasma recombination 
into highly excited atomic states. Treating the evolution of the plasma on the basis of kinetic 
equations, while ionization/excitation and recombination are incorporated using rate equations, we 
have investigated theoretically the Rydberg-to-plasma transition. Including the influence of spatial 
correlations on the plasma dynamics in an approximate way we find that ionic correlations change 
the results only quantitatively but not qualitatively. 
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Advances in cooling and trapping of neutral gases have 
opened up a new branch of atomic physics, namely dy- 
namics in ultracold (T C 1 K) systems. One interesting 
topic is the physics of ultracold Rydberg gases and ultra- 
cold neutral plasmas. In recent experiments, ultracold 
neutral plasmas have been produced from a small cloud 
of laser-cooled atoms confined in a magneto-optical trap 
01 0-JiJ' S 0- I n the experiments performed at NIST 
lii Li l| ! a plasma was produced by photoionizing laser- 
cooled Xe atoms with an initial ion temperature of about 
10 /iK. By tuning the frequency of the ionizing laser, the 
initial electron energy E e was varied corresponding to 
a temperature range IK < E e /k-Q < 1000 K, and the 
subsequent expansion of the plasma into the surround- 
ing vacuum was studied systematically. Remarkably, a 
significant amount of recombination was observed as the 
plasma expands, leading to the formation of Rydberg 
atoms from the plasma. In a complementary type of 
experiment 0, , ultracold atoms were laser-excited into 
high Rydberg states rather than directly ionized. In these 
experiments, also the reverse process has been observed, 
namely the spontaneous evolution of the Rydberg gas 
into a plasma. 

An issue in the theoretical considerations 0, 
stimulated by the experiments has been the question 
whether the evolving plasma is strongly coupled or not. 
Due to the very low initial temperature, the electron 
Coulomb coupling parameter T e (t — 0) is found to be 
significantly larger than one (T e = e 2 /(ak-QT e ), where 
a is the Wigner-Seitz radius). On the other hand, it 
has been pointed out in |(| that, for the initial condi- 
tions of the NIST experiments [3 0- Q. the develop- 
ment of equilibrium electron-electron correlations leads 
to a rapid heating of the electron gas, which brings T e 
down to order unity on a very short timescale. In Q, 
this was demonstrated by a molecular-dynamics simu- 
lation of the electron and ion motion. The calculation 
was limited to a short time interval (~ 100 nanoseconds) 
in the initial stage of the plasma evolution due to the 
large numerical effort required. For the first quantita- 
tive comparison with experiment the plasma dynamics 



has been modeled within a kinetic approach, while ion- 
ization, excitation and recombination has been treated 
by a separate set of rate equations 7]. The electron 
dynamics was described in an adiabatic approximation 
with the atom and ion temperatures set to zero. Since 
this model does only account for the mean-field poten- 
tial created by the charges, possible correlation heating 
could not be described. Nevertheless, due to heating by 
three-body recombination events, T e does not exceed a 
value of w 0.2 during the plasma expansion. Hence, the 
influence of electronic correlations on the plasma dynam- 
ics could be neglected. With zero ionic temperature the 
ionic Coulomb coupling parameter, however, is infinite 
in the framework of this model. Consequently, the role 
of ion-ion correlations could not be estimated. 

Our description is similar to that of Q with small 
differences, e.g., inclusion of black-body radiation as a 
source for photoionization. This is necessary to describe 
the initial ionization of the Rydberg gas in the experi- 
ments 0, Q since the ionization rate through cold atom- 
atom collisions is much lower at these densities. More 
importantly, in addition to Q we allow for possible spa- 
tial correlation effects. Briefly, we use a set of kinetic 
equations for the plasma expansion, combined with rate 
equations for the description of ionization/excitation and 
recombination. The kinetic equations are derived from 
the first equation of the BBGKY hierarchy, which yields 
the standard Vlasov equation for the evolution of a colli- 
sionless plasma, augmented by additional terms account- 
ing for spatial correlations. Since the electronic Coulomb 
coupling parameter is well below unity under all exper- 
imental conditions realized so far (as we have convinced 
ourselves including electronic correlations on the basis of 
two-component Debye-Hiickel theory; see also |7j) we ne- 
glect them and treat the ions as a one-component plasma 
embedded in a neutralizing background. 

Assuming that the correlation function vanishes at dis- 
tances larger than the correlation length a c and that the 
spatial density varies slowly over a distance a c , we find 
that the force induced by ionic correlations can be ex- 
pressed in terms of the local correlation energy per par- 
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where the spatial density pi is related to the one-particle 
distribution function /j by pi — J dvfi and g(y) is the 
correlation function. The correlation energy for a one- 
component plasma in thermal equilibrium is a function 
of the Coulomb coupling parameter only , and simple 
analytical formulae have been given for a wide range of 
T. We use 
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derived from the ion sphere model which is a good ap- 
proximation for r > 1 |ll| . 

The timescale for relaxation towards the individual 
equilibrium of electrons and ions, respectively, is deter- 
mined by the corresponding inverse plasma frequency 
Wp 1 = yj m /Airpe 2 ha. Under typical conditions of the 
experiments 0, Q> @>Wi ^ ^ s of the order of 10~ 9 s for elec- 
trons, while equilibration of the ions takes between 10~ 7 s 
and 10 _5 s depending on p. These timescales for equili- 
bration have to be seen relative to the expansion time of 
the plasma, which is of the order of several 10 _6 s. Hence, 
an adiabatic approximation for the electron distribution 
function can be safely applied, such that the kinetic equa- 
tion for the electronic component of the plasma yields a 
relation allowing one to express the mean-field potential 
in terms of the density p e , which can be inserted into the 
corresponding equation for the ions. 

The much slower equilibration of the ions renders 
an adiabatic approximation a priori difficult. However, 
even in situations where static parameters, e.g., relevant 
masses, periods or rates, speak against an adiabatic treat- 
ment it is sometimes made possible through dynamical 
adiabaticity (see, e.g., [HI). Of course, this can only 
be verified a posteriori by more elaborate calculations 
without adiabatic approximation. For now, we simply 
assume an adiabatic time evolution of the plasma, i.e. we 
estimate the ionic correlation energy from its equilibrium 
value eq. J5J) and use a homogeneous ion temperature at 
all times. This is the opposite limit to a zero ionic tem- 
perature which would never lead to equilibration. Hence, 
our results are expected to clarify whether (and to what 
extent) correlation effects can have any influence on the 
dynamics at all. 

The large difference in timescales between ionic and 
electronic motion justifies the quasineutral approxima- 
tion p e = pi = p 14] applicable under the experimental 
conditions in nj, |2j, y, |j] . Under these circumstances a 
closed equation for the ion distribution function is ob- 
tained and the resulting equation without the correlation 
term permits selfsimilar analytical solutions |m| . One of 
them is the gaussian profile describing the initial state in 



the experiments under consideration. Exact selfsimilar 
solutions exist even in the case of an additional linear 
force ^ij. This condition is satisfied by the correlation 
force to a good approximation save for the outer periph- 
ery of the plasma 16] . Hence, we can write 
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where a is the width of the spatial density distribution 
and w(r) = 7r is the hydrodynamic velocity (see also0)- 
Substitution of Eq. © finally leads to the set of equations 
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where W c = 1/3/ drpd(uup)/dp arises from the corre- 
lation pressure and M is the total mass of the plasma. 
The thermal energy is determined from the total energy 
of the system 
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i.e. the sum of the kinetic energy and the correlation en- 
ergy U c = J druup- The set of equations (QJ reduces to 
the one used in [7j if the correlation term is dropped and 
Tt = 0. More recently, it has been shown by compari- 
son with more sophisticated calculations that this simple 
ansatz provides a good description of the plasma dynam- 
ics even if some of the assumptions, such as the quasineu- 
trality of the plasma and the gaussian distribution of the 
densities, only hold approximately ^3 • 

The low-temperature enhancement of the expansion 
velocity of the ultracold plasma produced in can be 
well described by combining a hydrodynamic descrip- 
tion of the expansion dynamics with conventional rate 
equations accounting for inelastic collisions between the 
plasma particles and Rydberg atoms 7]. As a result, the 
total kinetic energy becomes a function of time since the 
Rydberg atoms act as energy sinks and sources. In our 
calculations, we include bound-bound transitions by elec- 
tron impact excitation and deexcitation, and bound-free 
transitions by three-body recombination, electron impact 
ionization, and black-body radiation. We use the col- 
lision rate coefficients derived by Mansbach and Keck 
[lif and black-body photoionization is described by first- 
order perturbation theory 0], using an as ymp totic ex- 
pression for the atomic oscillator strengths [20j. 

Finally, the question of the expansion velocity of the 
Rydberg atoms has to be addressed, since they are not 
driven by the Coulomb drag of the expanding electrons. 
Under typical experimental conditions (p e — 10 9 cm~ 3 , 
T e = 20 K), the timescale for electron-atom-collisions is 
of the order of 10 _8 s. Therefore, in the course of the ex- 
pansion of the plasma (10 _5 s) any given atom with bind- 
ing energy of the order of several k B T will recombine and 
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FIG. 1: Time evolution of the degree of ionization for the 
following initial conditions: density of atoms p = 8- 10 9 cm -3 , 
atomic temperature T a = 140 pK, plasma width a — 60 /im 
and initial principal quantum number of the Rydberg atoms 
no — 70. Results are shown with ionic correlations (solid) and 
without (dashed). 

re-ionize many times, changing its character between ion 
and Rydberg atom. In addition, collisions between ions 
and atoms significantly equilibrate the hydrodynamical 
velocities even for lower Rydberg states |2l|. Following 
this reasoning we assume equal hydrodynamical veloci- 
ties and density profiles for the ions and atoms. This 
implies that the expansion of the neutral Rydberg atoms 
can be simply taken into account by changing the mass 
of the ions to an effective mass M in Eq. (J2J which is the 
total mass of the total system (plasma + atoms). 

We have simulated the expansion of an initially fully 
ionized plasma as reported in 0, ■ In accordance with 
previous calculations 0, we find quantitative agreement 
for the electron energy dependence of the asymptotic ex- 
pansion velocity. Furthermore, the calculations qualita- 
tively reproduce the nonmonotonic time dependence of 
the number of detected atoms observed in |jg. We note 
that the ionic correlations do not significantly change the 
plasma dynamics for these experiments. 

The results of our calculations for the experiments 
0,0 are summarized in figures 1 to 3. Figure 1 shows the 
degree of ionization of the Rydberg gas, i.e. the plasma 
fraction of the system, as a function of time during the 
expansion of the system. For comparison, the result with- 
out correlation (U c = = in eqs. 0} and ©) is also 
shown. As has been observed experimentally, it takes of 
the order of two to three microseconds before significant 
plasma formation occurs which sets in with an ionization 
avalanche. Note that the early time development includ- 
ing the ionization avalanche is obviously not influenced 
by the ionic correlations. Their development requires a 
significant amount of charged plasma particles which are 
only available after the ionization avalanche. The pro- 
duction of initial "seed charges" by black-body radiation 
does not depend on ionic correlations. The ionization 
avalanche appears somewhat (w l^ts) earlier than in the 
experiment since in reality the first electrons produced 
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FIG. 2: Degree of ionization as a function of time for various 
initial conditions (no = 70). a) a = 80 /im fixed, p = 2.5 • 10 9 
(dotted), 5 • 10 9 (dashed), 1 ■ 10 10 cm" 3 (solid); b) p = 5 ■ 10 9 
cm" 3 fixed, a = 60 (solid), 80 (dashed), 100 pm (dotted); 
c) N a = 114000 fixed, a = 60 (solid), 80 (dashed), 100 pm 
(dotted). 

by the black-body radiation will leave the atom cloud 
until a sufficiently strong positive space charge builds up 
to trap the electrons. With the increasing number of free 
charges, the plasma forms quickly, followed by a partial 
back evolution into a Rydberg gas. Fig. shows that 
this part of the gas dynamics may indeed be affected by 
ion-ion correlation, which leads to a more efficient recom- 
bination in the final stage of the plasma evolution. The 
reason for this behavior can be understood from eqs. (0} 
and (0. Together, W c and U c lead to an additional ac- 
celeration (2/9U c if Eq. (J2J) is used for u) in addition 
to the ideal thermal electron pressure k^T e . Hence, the 
adiabatic electron cooling during the expansion is faster 
than without correlations. In turn recombination is ac- 
celerated, since the corresponding rate strongly increases 
with decreasing temperature. We note that correspond- 
ing effects of the same order of magnitude are also found 
for the experiments 0, 0] • The reason that they do not 
significantly affect the expansion dynamics of the plasma 
part lies in the fact that the correlations mainly influence 
the recombination into high-lying states at later stages of 
the evolution. Since these are very weakly bound, they 
hardly influence the kinetic energy of the system which 
determines the expansion velocity. Hence, good agree- 
ment is found between calculations with and without in- 
clusion of correlation effects as long as quantities related 
to the plasma expansion are compared. 

Figure [21 illustrates the dependence of the plasma evo- 
lution on the initial conditions. For fixed initial width 
of the Rydberg atom cloud (Fig. |2Ji), the plasma for- 
mation occurs earlier with increasing density, due to the 
fact that ionization is faster and more efficient at higher 
densities. The later stages of the evolution, on the other 
hand, depend only weakly on the density. This may be 
attributed to the fact that in the evolution equation 0}, 
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FIG. 3: Level distribution of Rydberg atoms after lps, 1.4/zs, 
4/is and 25/LtS. 



lution. Microscopic molecular-dynamics simulations are 
limited to much shorter timescales due to the large nu- 
merical effort required. We have estimated the maximum 
effect of ionic correlations by including them into the 
model assuming instant equilibration of the ions, fn this 
respect the calculations without equilibration (no corre- 
lation, ionic temperature zero) and with instant equili- 
bration should bracket the actual ionization yield of the 
Rydberg gas. Given the difference in the ionization yields 
with and without ionic correlation in fig. ^ more elabo- 
rate calculations are worthwhile although at present it is 
unclear whether the plasma evolution including correla- 
tion can be followed over sufficently long times, e.g., with 
P 3 M-codes. 
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p enters only in the form N/M, i.e. essentially as the 
degree of ionization. For fixed initial density p, the ion- 
ization avalanche occurs at the same time, in accordance 
with the argument made previously (Fig. 03). However, 
for smaller width, the plasma expands faster, leading ul- 
timately to lower electron temperatures and increased re- 
combination. Finally, if the number of atoms N a is kept 
fixed the density changes when the size of the atom cloud 
is varied. Hence, Fig. |2t can be interpreted as the com- 
bined effect seen in Figs. ^.,b: The ionization avalanche 
occurs at different times for each curve (due to different 
atomic densities) and the relaxation leads to different fi- 
nal degrees of ionization (due to the different size of the 
atomic cloud). 

Figure |21 shows the non-ionized part of the gas, more 
precisely the population of Rydberg levels in the initial 
and final stage of the expansion, as well as at the time 
of the maximum degree of ionization. Initially, all atoms 
are prepared with a principal quantum number uq = 70. 
At later times, the interplay between ionization, recom- 
bination and exciting and deexciting collisions leads to 
a decrease of the average excitation. This decrease in 
excitation is the main source of energy that triggers the 
ionization and the plasma expansion, energy absorbed 
from the black-body radiation field is negligible in com- 
parison. Finally, recombination repopulates higher levels 
(which "freeze out" due to the rapid decrease of elec- 
tron temperature), but the peak of the distribution is 
still found at relatively low excitations. It would be very 
interesting to measure the final distribution of Rydberg 
atoms, which should provide a stringent test for theoret- 
ical models. 

In summary, we have described the formation of a cold 
plasma from a highly excited Rydberg gas, using a sim- 
ple model based on kinetic equations for the plasma evo- 
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